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Abstract 

The Solar Chimney Power Plant System (SCPPS) is a simple 
solar thermal power plant that is capable of converting solar 
energy into thermal energy in the solar collector. In the 
second stage, the generated thermal energy is converted into 
kinetic energy in the chimney and ultimately into electric 
energy using a combination of a wind turbine and a 
generator. This study is to conduct a more detailed 
numerical analysis of solar chimney power plant system 
with a curved junction. This paper summarizes a numerical 
study flow in a solar chimney with a curved junction. The 
fluid is the air (Pr=0.702), and considered as a Newtonian 
and incompressible fluid. The governing equations are taken 
to be in the vorticity-stream function formulation in 
hyperbolic coordinates by using the Boussinesq 
approximation. For heating conditions, an isothermal walls 
of the collector has been supposed. Solution of the defined 
equations has been done with numerical control volume 
method. Additionally, the effect of the system geometry (the 
curved junction) on the heat transfer phenomenon in the 
solar chimney has been examined. Results are related to the 
temperature distribution and the velocity field in the 
chimney and in the collector. 
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Introduction 

A Solar Chimney Power Plant System (SCPPS) is 
composed of a solar collector, to raise the energy level 
of the air by greenhouse effect of a chimney tower 
ensuring the circulation of air per gradient of density, 
and of an aero generator to produce electric power. 

Professor Schlaich of Stuttgart originally proposed the 
solar chimney concept in the late 1970s (Schlaich, 1995). 
Less than 4 years after he presented this ides at a 
conference, construction on a pilot plant began in 
Manzanares, Spain, as a result of a joint venture 
between the German government and a Spanish utility. 
A 36 kW pilot plant was built, which produced 



electricity for 7 years, thus proving the efficiency and 
reliability of this novel technology. The chimney tower 
was 194.6m height, and the collector had a radius of 
122m. Fundamental investigations for the Spanish 
system were reported by (Haaf, et al., 1983) in which 
a brief discussion of the energy balance, design criteria, 
and cost analysis was presented. Efforts were focused 
essentially on analyzing performances and cost of 
solar chimney. (Bernardes, et al., 1999) presented a 
theoretical analysis of a solar chimney, operating on 
natural laminar convection in steady state. In order to 
predict thermo-hydrodynamic behavior of air, 
temperature conditions were imposed on entrance, so 
as to guarantee steady laminar flow along the device. 
The mathematical model was analyzed by the method 
of Finite volumes in generalized coordinates. Velocity 
field and temperature distribution in the flow were 
obtained under imposed thermal conditions. (Von 
Backstrom & Fluri, 2006) investigated analytically the 
validity and applicability of the assumption that, for 
maximum fluid power, the optimum ratio of turbine 
pressure drop to pressure potential is 2/3. (Von 
Backstrom & Gannon, 2000) were interested mainly 
in a one-dimensional compressible flow for the 
thermodynamic variable as dependence on chimney 
height, wall friction, additional losses, internal drag 
and area exchange. (Pretorius & Kroger, 2006) 
evaluated the influence of a developed convective heat 
transfer equation, more accurate turbine inlet loss 
coefficient, quality collector roof glass and various 
types of soil on the performance of a large scale solar 
chimney power plant. (Ming, et al., 2006) presented a 
mathematical model to evaluate the relative static 
pressure and driving force of the solar chimney power 
plant system and verified the model with numerical 
simulations. (Maia et al.,2009) presented a theoretical 
analysis of a turbulent flow inside a solar chimney. 
They showed that the most important physical 
elements in a solar chimney system are the tower 
dimensions as they cause the most significant 
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variation in the flow behavior. An increase in the 
height and diameter of the tower produces the rise in 
the mass flow rate and decrement in the flow 
temperature. (Koonsrisuk & Chitsomboon, 2009) 
studied about the ventilation efficiency of solar 
chimney by comparing between five of the 
mathematical simulations and five of CFD simulation 
both from the previous researches. (Koonsrisuk & 
Chitsomboon, 2007) proposed dimensionless variables 
to guide the experimental study of flow in a small- 
scale solar chimney: a solar power plant to generate 
electricity. Computational fluid dynamics (CFD) 
methodology was employed to obtain results that are 
used to prove the similarity of the proposed 
dimensionless variables. (Pastohr et al., 2004) 
presented a numerical simulation result in which the 
storage layer was regarded as solid. In their paper, 
conjugate numerical simulations of the energy storage 
layer, the collector and the chimney have been 
conducted, and the characteristics of the heat storage 
system, and the flow and heat transfer in the whole 
system have been studied. (Zhou, et al., 2007) have 
performed an experimental study in a solar chimney. 
A pilot experimental solar chimney power setup 
consisting of an air collector of 10 m in diameter and 
an 8 m height chimney was built. The authors noted 
that the temperature difference between the collector 
outlet temperature and the one of the ambient usually 
might reach as much as 24.1 8 °C, which generates the 
driving force of the air flow in the setup. Their data 
analysis showed an air temperature inversion in the 
latter chimney after sunrise and this is due to the 
increase of solar radiation from the minimum. The 
phenomenon disappears once a driving force is 
generated by a temperature high enough to overcome 
it. (Koonsrisuk et al., 2010) described the constructal- 
theory search for the geometry of a solar chimney. 
(Sangi et al., 2011) developed a mathematical model to 
accurately describe the solar chimney power plant 
mechanism. Numerical profiles for the temperature, 
velocity and pressure in the collector of the solar 
chimney power plant were shown for three different 
solar radiations. (Chergui et al., 2010) presented a 
study considering the heat transfer process and the 
fluid flow in the collector and the chimney under some 
imposed operational conditions. The temperature 
difference through the Rayleigh number, the velocity 
field and the temperature distribution, through the 
system, have been evaluated. Results showed local 
flow characteristics and it appears that for most 
Rayleigh numbers, the flow seems laminar except for 



Raleigh number of 10 8 where there are some 
disturbances. 

Problem Formulation and Basic Equations 

A solar chimney has been taken into consideration 
which contains an incompressible fluid (air). The roof 
and ground surface create a vertical temperature 
gradient (active walls) (Fig.l). 

The physical properties of the fluid are constant, apart 
from the density q whose variations are at the origin of 
the natural convection. Viscous dissipation is 
neglected, just as the radiation (emissive properties of 
the two walls being neglected). It was accepted that 
the problem is bidimensionnal, permanent and 
laminar. 



£1" dv £u 




FIG. 1 PROBLEM STUDY AND BOUNDARY CONDITIONS 

The laminar natural convection equations within the 
framework of the Boussinesq approximation are 
written: 



Continuity equation 



du dv „ 
— + — = 

dx dy 



Vorticity equation 
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dx 

Energy equation 
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(3) 



It is convenient to define a reference frame such as the 
limits of the system result in constant values of the 
coordinates. The passage of the Cartesian coordinates 
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(x, y) to the hyperbolic coordinates (rj, 6>) is obtained 
by the following relations: 



r + r/ 
2 



y = 



r-tj 



(4) 



With: 



r = ^ 2 +0 2 (5) 
Equations (1), (2) and (3) written then respectively: 

(6) 



„ 57* T 8T I A 

V„ +Vn = 

n drj 9 30 h pCp 



V, dco V„ dco _ of 3 2 co 5 J » 
h 3r/ + h d8~ h\dr/ 1+ d6 2 



( ' *2 



2r\ 



3 A T d z T 
drj 2 30 2 



(7) 



r + ^ar + r-77 (S) 
2r dt] V 2r 50 I 



Where, /i is the metric coefficient. 
The equation for the stream function: 



a 2 [e^ 2 56> 2 



(9) 



The quantities characteristic used for the 
dimensionless problem are the AT = Ti, - Tc between 
The roof and ground surface, the metric coefficients in 
hyperbolic coordinates (h) as reference length and the 
thermal diffusivity of fluid a as characteristic velocity. 
The dimensionless mathematical model obtained is: 



1 
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The boundary conditions are the following ones: 
The roof of the collector 

V* = V g *= ^ = ^ = 0andr = 



(10) 



(11) 



(12) 



(13) 



The axis of symmetry 



38 



a 2g. = 0i v e - = >£ = o,°£ = o 

d8 2 u drj d8 



The ground 



V* = V e *= 55! = 2*1 = and r = 1 

^ w 58 drj 



(14) 



(15) 



(16) 



Collector inlet 



= "e = -rr = ~r~ = and T 
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= 



Chimney outlet 
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= (18) 



Numerical Method 



To solve this system of equations with associated 
boundary conditions equations, a numerical solution 
by the method of finite volumes has been taken into 
account, presented by (Patankar, 1980) using a 
numerical code (computer programme). The power 
law scheme was used for the discretization. The 
iterative method used for the numerical solution of 
algebraic system of equations (Matrix) is the Gauss- 
Seidel, with an under-relaxation process. The 
following condition is used to declare convergence: 

|maxr n+1 -maxr n | 



< 10' 



(19) 



I maxT n+1 

Where: n+1 is the current iteration and n is the 
previous iteration. 

Fig. 2 shows the physical domain and the computational 
domain. 




The ajds of symmetry 



Tile roof of tile collector 



%. The roof of the collector 





^ 1 The axis of symmetry The ground 



The ground 

FIG. 2 PHYSICAL DOMAIN AND COMPUTATIONAL DOMAIN 

Results and Discussion 

The grid dependence has been investigated using 
different mesh sizes before settling to a mesh size of 
(350 x 22). 

Our objective is to analyze the effect of geometry on 
heat transfer and flow of air into the chimney. For this 
reason, we presented the isotherms and isovelocity 
lines in two different geometries for three different 
Rayleigh numbers. 

For the validation of the computational problem, our 
results were compared with those of literature: (Sangi 
et al., 2011), and the study of (Pastohr et al., 2004) and 
(Chergui et al., 2010). 

Figs. 3-5 show the dimensionless isothermal lines for 
the first geometry and for Rayleigh numbers equal to 
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10 3 , 10 4 and 10 5 respectively. When the Rayleigh 
number is small, as being lower or equal to 10 4 , the 
heat transfer is essentially conductive, so the isotherms 
have the same form as the walls and the maximum 
temperature is located near the ground in the collector 
due to heat transfer exchange between this surface and 
the airflow beneath the cover. For a Rayleigh number 
equal to 10 s swirls in isothermal lines are also 
observed in the collector area. 




FIG. 3 ISOTHERMAL LINES FOR THE FIRST GEOMETRY AND 
Ra=10 3 




FIG. 4 ISOTHERMAL LINES FOR THE FIRST GEOMETRY AND 
Ra=10 4 




FIG. 5 ISOTHERMAL LINES FOR THE FIRST GEOMETRY AND 
Ra=10 5 




FIG. 6 ISOVELOCITY LINES FOR THE FIRST GEOMETRY AND 
Ra=10 3 




FIG. 7 ISOVELOCITY LINES FOR THE FIRST GEOMETRY AND 
Ra=10 4 




FIG. 8 ISOVELOCITY LINES FOR THE FIRST GEOMETRY AND 
Ra=10 5 

Figs. 6-8 represent the isovelocity lines for the first 
geometry and different Rayleigh number values. For 
Ra = 10 3 and 10 4 , the velocity magnitude increases with 
the rise of the Rayleigh number and its maximum 
value is located approximately at the inlet of the 
chimney as it is reported in literature (Pastohr et al., 
2004), (Sangi et al., 2011) and (Chergui et al., 2010). For 
these Rayleigh numbers, isovelocity lines are smooth 
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characterizing a laminar flow, along the chimney and 
the collector, except for the region of curved junction 
and the inlet of the collector where different zones of 
flow recirculation can be observed due to instability of 
airflow in these areas. For Ra = 10 s an instability 
(different zones of flow recirculation) of airflow along 
the collector can be observed. The isovelocity lines 
values show an appreciable increase in the flow. This 
means that the transfer is done primarily by 
convection and predominates on the conduction. 

When the distance between the ground and the roof is 
increased (second geometry), a disturbance of 
isothermal lines can be observed in the collector zone 
at the beginning of a Rayleigh number equal to 10 4 
(Figs. 9-11) and the natural convection is dominant. 

Figs. 12-14 illustrate the dimensionless velocity lines 
for second geometry and for Rayleigh 10 3 ,10 4 and 10 s 
respectively. It should be noted that its value is higher 
than the one calculated in the first geometry and 
always its maximum is located at the inlet of the 
chimney. A large zone of flow recirculation is inverted 
located below the curved wall. 




FIG. 9 ISOTHERMAL LINES FOR THE SECOND GEOMETRY AND 

Ra=10 3 




FIG. 10 ISOTHERMAL LINES FOR THE FIRST GEOMETRY AND 
Ra=10 4 




FIG. 11 ISOTHERMAL LINES FOR THE FIRST GEOMETRY AND 
Ra=10 5 




FIG. 12 ISOVELOCITY LINES FOR THE SECOND GEOMETRY 
AND Ra=10 3 




FIG. 13 ISOVELOCITY LINES FOR THE SECOND GEOMETRY 
AND Ra=10 4 




FIG. 14 ISOVELOCITY LINES FOR THE SECOND GEOMETRY 
AND Ra=10 5 
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Fig. 15 shows the air velocity profile through the 
collector for the two geometries considered and 
Rayleigh 10 3 . The velocity increases through the 
collector by decreasing the radius, but it increases 
more sharply by reaching the chimney base. When the 
collector radius is constant, an increase of distance 
between the ground and the roof (second geometry) 
causes an increase of the air velocity. The same results 
are observed in literatures by (Pastohr et al., 2004), 
(Sangi et al, 2011) and (Chergui et al., 2010). 
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FIG. 15 VELOCITY PROFILE THROUGH THE COLLECTOR FOR 
THE TWO GEOMETRIES CONSIDERED AND Ra=10 3 

Conclusions 

In this paper a thermo-hydrodynamic analysis for air 
motion in natural convection, laminar flow and steady 
state has been presented for a solar chimney with 
prescribed boundary conditions. A validated computer 
program was adapted to the solar chimney configuration 
to solve the governing equations, using the method of 
finite volumes. Meanwhile, the effect of the system 
geometry with a curved junction on the heat transfer 
process and the fluid flow in the system has been 
investigated, along with two different geometries 
considered. The simulations were executed based on 
three values of the Rayleigh number: Ra=10 3 , Ra=10 4 and 
Ra=10 5 . 

Results showed that the maximum velocities are 
gotten at the inlet of the chimney tower and and its 
values were increased by increasing the difference 
between ground and roof of the collector for all 
Rayleigh numbers. 

This result will help the solar chimney designer 
correctly locate the turbine in the solar chimney power 
plant. 

Nomenclature 

Symbols 

Cp: Specific heat at constant pressure, J kg- 1 K 1 ; 
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h: 


Metric coefficient, m; 


H: 


Dimensionless metric coefficient; 


P<: 


Prandtl number; 


Ra: 




1 . 


T-Hiiirl'c t~£*mn£*rz)t~i ir£* Tv~* 
1 1UH-1 & LCllllJCl CllUl c, XV, 


Tv. 


Temperature of the ground, K; 


Tc. 


Temperature of the roof, K; 


AT: 


Th- Tc, K; 


u , v: 


Velocities components according to coordinates 



x and y, m s _1 ; 



V n , V e : Velocities components according to coordinates 
r\ and 8, m s 1 ; 

x, y: Cartesian coordinates, m; 
Greek Symbols 





Thermal expansion coefficient, K 1 ; 


X: 


Thermal conductivity, W.rrr 1 K 1 ; 


u: 


Kinematic viscosity, m 2 s 1 ; 


P- 


Density, kg.nr 3 ; 


11,0: 


Hyperbolic coordinates; 


y: 


Stream function, m 2 s 1 ; 


to: 


Vorticity ,s _1 ; 


Subscripts 


h : 


Hot; 


c: 


Cold; 




Dimensionless parameters. 
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